Invariant number triangles, eigentriangles and 

Sonios-4 sequences 



Paul Barry 

School of Science 

Waterford Institute of Technology 

Ireland 

^ ; pbarrySwit . ie 

^ Abstract 

t^ ; 

^N , Using the language of Riordan arrays, we look at two related iterative processes on 

matrices and determine which matrices are invariant under these processes. In a special 

^j ' case, the invariant sequences that arise are conjectured to have Hankel transforms that 

\^ • obey Somos-4 recurrences. A notion of eigentriangle for a number triangle emerges and 

r~| ! examples are given, including a construction of the Takeuchi numbers. 

1 Introduction 



o 



% 



^ . In this note, we shall define transformations on invertible lower-triangular matrices involv- 

O ! ing the down-shifting of elements and taking an inverse. The invariant matrices for these 

transformations turn out to be simple Riordan arrays [9], with generating functions easily 
lO I described by continued fractions [4, 13]. These matrices have close links to the Catalan 

t^^ ■ numbers C„ = ^^ in)- ^^ ^^^ ^^^^ °^ ^ particular two-parameter transformation, special 

sequences defined by this process appear to have Hankel transforms that satisfy Somos-4 
type recurrences [3]. Again using Riordan arrays we can characterize these sequences. 

We recall that the Riordan group [9, 11], is a set of infinite lower-triangular integer 
matrices, where each matrix is defined by a pair of generating functions g{x) = 1 + gix + 
g2x'^ + . . . and /(x) = fiX + f2X^ + . . . where /i ^ [11]. The associated matrix is the 
matrix whose i-th column is generated by g{x)f{xy (the first column being indexed by 0). 
The matrix corresponding to the pair f,g is denoted by ((?, /). The group law is then given 

by 

{g,f)-{h,l) = {g{hof),lof). 

The identity for this law is / = (l,x) and the inverse of {g,f) is {g,f)^^ = (l/(fi' o /),/) 
where / is the compositional inverse of /. This is also called the (series) reversion of /. A 
Riordan array of the form {g{x),x), where g{x) is the generating function of the sequence a„, 
is called the sequence array of the sequence a„. Its general term is a„_fc, or more accurately 
[k < n]an-k (where [P] is the Iverson bracket [6], defined by [P] = 1 if the proposition V is 
true, and [P] = if P is false). Such arrays are also called Appell arrays as they form the 
elements of the so called Appell subgroup. 



If M is the matrix {g, /), and a = (oq, ai, . . .)' is an integer sequence with ordinary gener- 
ating function A (x), then the sequence Ma has ordinary generating function g{x)A{f{x)). 
The (infinite) matrix [g, f) can thus be considered to act on the ring of integer sequences 
Z^ by multiphcation, where a sequence is regarded as a (infinite) column vector. We can 
extend this action to the ring of power series Z[[x]] by 

{g, /) : A{x) -^ {g, f) ■ A{x) = g{x)A{f{x)). 

Example 1. The binomial matrix B is the element (y3^) yz^) of the Riordan group. It 
has general element (^). More generally, B™ is the element (jz^i j^^) of the Riordan 
group, with general term (2^m"~^. It is easy to show that the inverse B~"^ of B™ is given 

by (t-^,tt^). 

In the sequel, we shall assume that all matrices and sequences are integer valued. 



2 The (a, 6)-Process 



We start by defining an operation on lower-triangular matrices which have I's on the diag- 
onal. Thus let M be of the form 
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Then we take the inverse M{a, b) ^ of this matrix. Let us call this process the (a, 6)-process. 
We have the following proposition. 



Proposition 2. Let f{x) be the power series defined by 

1 



/(^) 



1 — ax — {b — l)x^ — x^/(x) 



(3) 



Then the Riordan array 

is invariant under the [a, b)- operation. 

Proof. By equation (3), we see that f{x) = Yl'ilo^^i^^ where ao = 1. Then 

x^ f{x) = x'^ao + x^ y^ cLi+ix^ = x"^ + x^ y^ ca+ix^- 

i=0 j=0 

We obtain 

l — ax — {h — l)x^ — x'^f{x) = l — ax — bx'^ + x^ — x'^—x^y^ai+ix^ = 1 — ax — bx'^ — x^y ^ai^ix\ 

i=0 i=0 

Thus we wish to prove that 

{f{x),x) = (1 — ax — {b — l)x^ — x'^ f (x) , x)~^ , 
or equivalently that 

{f{x),x)~^ = (1 — ax — {b — l)x^ — x^/(x), x). 
Now 



and hence we wish to estabhsh that 

1 



1 — ax — {b — l)x — X f{x). 



fix) 

But this follows immediately from the definition of /. D 

Let a„ denote the n-th element of the first column of (/(x), x). Then the {n, k)-th element 
of (/(x),x) is given by 

[k < n]an-k- 

Thus we need only a knowledge of a^ to describe all elements of the matrix. 

Proposition 3. Let 

9{x) 



1 — ax — bx'^ — x'^g{x) 
Then 

where C„ = -r-^i ") is the n-th Catalan number A000108. 



Proof. Solving the equation 



1 — ax — hx'^ — x'^g{x) 
gives us 



1-ax-bx'^ - Jl - 2ax + (a^ - 26 - 4)x^ + 2abx^ + b^x^ 
9{x) = ga,b{^) = ^2 

With this value, we then have the Riordan array factorization 

/ /NX f 1 X \ ( ( 2n 9aA^) 



1 — ax — bx"^ ' 1 — ax — bx 



X 



1 ^^/'l ^^/^/2^ 9a,b{x) 



^ 1 — bx"^ ' 1 — bx'^ J \1 — ax' 1 — ax / \ ' x 

where 

, , 1 - VI - 4x 

c[x) = 

^ ' 2x 

is the g.f. of the Catalan numbers, and c(x^) is the g.f. of the aerated Catalan numbers 
1,0,1,0,2,0,5,0,.... Thus 

1 X \ f 1 X \ , 2\ 



D 



The result follows from this. 
Corollary 4. 

A;=0 ^ ^ i=0 \ -^ / 

We note that if we start with any matrix of the form (1), and iterate the (a, 6)-process 
on it, then the limit matrix is {f{x),x). Thus the element of the Appell subgroup of the 
Riordan group (/(x),x) where 

1 
f{x) = ^ , 

X 

1 — ax — {b — l)x^ — ■ 



1 — ax — (6 — l)x'^ 



1 
is a "universal element" for the (a, 6)-process. 



3 A Somos-4 conjecture 

We have the following Somos-4 conjecture. 

Conjecture 5. The Hankel transform of the sequence an is a (a^,fe^ — a^) Somos-A sequence. 
By this we mean that the sequence hn of Hankel determinants 

'''n \(^i+j\0<i,j<n 

satisfies an (a, (3) Somos-4 relation 

ahn^ihn-3 + I3hl_2 
hn = 7 , n > 3, 

where a = a"^ and (3 = h'^ — a^ . 

Equivalently the Hankel transform of the sequence with general term 

k=0 ^ ^ j=0 ^ -J ^ 

is (conjectured to be) a (a^, (6 + 1)^ — a^) Somos-4 sequence. 

Example 6. We let a = 6 = 1. Then a„ is the sequence A128720 

1,1,3, 6, 16, 40, 109, 297, 836, 2377, 6869 . . . 

which counts the number of skew Dyck paths of semi-length n with no UUU's. The Hankel 
transform of this sequence is the (1,3) Somos-4 sequence A174168 which begins 

1, 2, 5, 17, 109, 706, 9529, 149057, 3464585, 141172802, 5987285341, .... 

Example 7. We take a = 1, 6 = 2 to get the sequence A174171 which begins 

1, 1, 4, 8, 25, 65, 197, 571, 1753, 5351, 16746 . . . , 

with (1,8) Somos-4 Hankel transform 

1, 3, 11, 83, 1217, 22833, 1249441, 68570323, 11548470571, 2279343327171, .... 

This is A097495 , or the even-indexed terms of the Somos-5 sequence. 

Example 8. We let a = 2, and b = —1. Then a„ is the sequence A187256 which begins 

1, 2, 4, 10, 28, 82, 248, 770, 2440, 7858, 25644, .... 

This sequence counts peakless Motzkin paths where the level steps come in two colours 
(Deutsch). The Hankel transform of this sequence is the Somos-4 variant A162547 that 
begins 

1, 0, -4, -16, -64, 0, 4096, 65536, 1048576, 0, -1073741824, .... 



4 The "(a)-process" and Narayana numbers 

We now look at the simpler "(a)-process" , whereby we send the matrix 



M 



/ 1 

m2,i 10 

m^^i ?7i3^2 10 

7714^1 7714^2 ^^14^3 1 

m^,i m^,2 "^5,3 m^A 1 

"^6,1 "^6,2 "^6,3 "^6,4 "^6,5 1 



\ 



to the matrix 
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-"^4,3 
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-ms^i -m5,2 -"^5,3 -"^5,4 -a 1 



V 



/ 



and then take the inverse to obtain M^ ^. We have the following result. 
Proposition 9. Let f{x) be the power series defined by 

1 



fix) 



Then the Riordan array 



1 — (a — l)x — xf{x) 
{f{x),x) 



is invariant under the [a) -process. 
Proof. We wish to show that 

(/(x), x) = (1 - (a - l)x - xf{x), x)~\ 
or equivalently that 

1 



ifix),xr' 



fix) 



,x 



(1 — (a — l)x — x/, x). 



But this follows immediately since by definition 



fix) 



(a — l)x — xf{x) 



(4) 



D 



We now remark that the continued fraction 

fix) = 

1 — (a — l)x - 



X 



1 — (a — l)x — 



X 



is the generating function of the Narayana polynomials A/'„(a) = ^22=0 ^n,kof' [1, 2, 12] where 
the matrix {Nn^k) is the matrix of Narayana numbers A090181 



/ 1 

10 

Oil 00 

13 1 

16 6 10 

1 10 20 10 1 



v ■■ 



\ 



J 



Hence the terms of the first column of {f{x),x) are precisely the Narayana polynomials in 



a: 



K. 



A:=0 



,fca 



In particular, for a = 1, we get 



G„ 



the Catalan numbers. 

As before, we note that if we start from an arbitrary matrix of the form Eq. (4), and 
iterate the (a)-process, then the limit matrix is {f{x),x). In particular, if a = 1, the limit 
matrix is the Catalan numbers sequence array (C„_fc): 

/ 1 

1 10 

2 110 

5 2 110 

14 5 2 1 1 

42 14 5 2 1 1 



V 



\ 



/ 



This is the Riordan array (c(x),x). 
By solving the equation 



fix) 



we see that 



ifix),x) 



1 — {a — l)x — xf{x) 

1 - (a - l)x - v/l-2(a+l)x + (a-l)2x2 
2x 
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which by the above is the matrix with (n, k)-th term 

[k < n]A/'„__fe(a). 



5 Eigentriangles 

We also have the following result. 
Proposition 10. Let M be a matrix 
By this we mean that if 



as in 



M{^ 



Eq. (1). Then M^ ^ is an eigentriangle of M. 

\ 



I I 

1 1 

rs.i 110 

^"4,1 ^"4,2 1 10 

^"5,1 ^"5,2 ^^5,3 110 

^"6,1 ^"6,2 ^^6,3 ^^6,4 1 1 



then 



MM^^ = 



\ 



/ 1 

r3,i 1 

r-4,1 r4,2 10 

^5,1 ^^5,2 r-5,3 10 

^6,1 ^^6,2 ^6,3 ^6,4 1 

rj,! r-y^a ^^7,3 ^7,4 ^^7,5 1 



/ 








V 



Note that the first column of M^ ^ is then an eigensequence of M 

Proof. We have 

MM{^ = I 

and hence 

fc-i 



\ 



/ 



Then for k ^ I, we have 



fe-i 
- ^ mk^ijTj^i + rk,i = for A; 7^ /. 



fe-i 

j=0 



Thus the (fc - 1, /)-th element of MM{^ is r^,;. 



(5) 



D 



Example 1 1 . The eigentriangle of the binomial matrix ( (^) ) is given by 
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31 
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1 ... 



E 



J 

where the first column entries are the Bell numbers. We note in passing that the production 
matrix [5] of the matrix E is equal to 

/ 1 1 
1 10 
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1 





1 
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1 





1 





15 
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1 





1 


52 


31 


14 


5 


1 






\ 



In this case, we have 



or 



the Bell numbers AOOOllO. 



n-l 
fc=0 



n — 1 
k 



Ofc, n > 0, ao 



Bell{n) 



Example 12. The eigentriangle of the skew binomial matrix ((„_^.)) is given by 

/ 1 ... \ 

1 10 

1 110 

2 2 110 
4 4 3 110 
11 11 7 4 1 1 



E 



\ 



where the first column 



1, 1, 1, 2, 4, 11, 33, 114, 438, 1845, 8458, . . . 
or A127782 is thus an eigensequence of {{^i^) (remark by Gary W. Adamson). We have 

n-l 



fc=0 



n 



k 
k 



Ofc, n > 0, ao = 1. 



Example 13. The eigentriangle of the sequence array for the Motzkin numbers M„ (i.e., 

the matrix with (n, k)-th term [k < n]Mn-k where M„ = Ylik=o {^k)^^) ^^ ^^^ sequence array 
for the sequence A005773 of directed animals An of size n. Thus 



n-l 



An = ^Mn-k~lAk. 



k=0 



We can characterize the eigentriangle E = {E{n, k)) corresponding to a matrix A = {A{n, k)) 
as follows. We define 



n-l 



E{n,j) = Y,Mn-l+J,k + j)E{k,j), with E{0,j) = l. (6) 



fc=0 

Then 



E{n, k) = [k< n\E{n - k, k). 



6 The Takeuchi numbers 

The Takeuchi numbers t„ A000651 are an example of a sequence that can be defined with 
the aid of the eigentriangle of the Catalan triangle {c{x),xc{x)) A033184 . We let T(x) be 
the generating function of the Takeuchi numbers. Our point of departure is (4) in [8]: 

^/ N cix) — \ x{2 — c{x))^, , ,, 

nx) = ^ + ^^J' t{xc{x)). 

1 — X V 1 — 4x 

We now note that 

(2-cW) , , 

Tr^lT = '=<"'■ 

so that [8] (4) becomes 

T{x) = — h xc{x)T{xc{x)). 

A. Jy 

In terms of Riordan arrays, we may write this as 

q(x] — 1 

((l,x) — (xc(x),xc(x))).T(x) = . 

1 — X 

Now while the matrix 

(l,x) — {xc{x) , xc{x)) 

is not a Riordan array, it is a special type of invertible matrix. The theory of eigentriangles 
tells us that its inverse is the eigentriangle of the Catalan matrix 

{c{x) , xc{x)) . 
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This eigentriangle begins 
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We then have 



/ 1 

1 
2 

5 




1 
2 

5 









1 



1 

3 



14 14 9 4 10 
42 42 28 14 5 1 



/ 1 

1 
2 
6 




1 
1 
3 





1 
1 

22 11 4 
92 46 17 









1 
1 1 
5 1 1 



/ V '■■ 
( 1 

2 
6 

22 
92 





1 

3 

11 

46 






1 
4 
17 







1 

5 











1 



426 213 79 24 6 1 



v 



c{x)-\ 



The sequence with g.f. -3-3^ is the sequence A014138 with general term 



k=0 



k+l, 



and thus the Takeuchi numbers are the image of this sequence by E. Now in this case A of 
Eq. (6) is the matrix {c{x), xc{x)) with (n, A;)-th term 

'2n - k\ k + 1 



A{n,k) 



n — k J n + 1 



Thus we get 



n-l 



%w)=Er":"/'"')^i^^(*.j). w,th E(o,,)=i, 



k=0 



n — 1 — k 



n + j 



and so 



fc-i 



k=0 j=0 

We note that the first column of E is essentially A091768 . 
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